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Abstract 

We predict the existence of an energy barrier for collapse in a system of two tunnel-coupled 
repulsive and attractive quasi two-dimensional condensates trapped in a double-well potential. 
The ground state in such a system can have a lower energy than that of the collapsing state. We 
find such ground states numerically and analytically determine the domain of their existence, they 
have a much larger share of atoms confined in the repulsive condensate. The energy barrier for 
collapse is due to the fact that the energy of the system increases if atoms tunnel from the repulsive 
condensate to the attractive one. 
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I. INTRODUCTION 



Bose-Einstein condensates (BEC) of degenerate quantum gases at zero temperature are 
described with a good accuracy by the mean-field theory based on the nonhnear Schrodinger 
equation with an external potential, i.e., the well-known Gross-Pitaevskii (GP) equation l| 
for the order parameter. The appearance of the order parameter is a consequence of the 
macroscopic quantum coherence of BEC, which was demonstrated experimentally and 
explained theoretically [4] with the use of the GP equation (see also the review Ref. [5|). 

The interplay between the quantum coherence and nonlinearity in BEC is a reach source 
of interesting new phenomena. In this paper we predict such a new phenomenon related to 
the collapse instability of an attractive BEC. The collapse in the condensate with a negative 
scattering length is extensively studied experimentally and theoretically (see, for instance. 



Refs. 



SB 



IG, 



111 ] and the references therein). The multi-component (multi-species) 
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BECs were also shown to suffer from the collapse instability [12, ll3|, [IJ, ll5| . Collapse of 
a repulsive BEC can be induced by a sudden change of its positive scattering len gth to a 
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negative value by application of the magnetic field near the Feshbach resonance 
Here we show that with the use of a coherent coupling between two BECs with the opposite 
interactions (one repulsive and the other attractive) it is possible to modify significantly the 
properties of the collapse instability. An experimental realization is provided by a pair of 
condensates trapped in the double-well potential with far separated wells, when the atomic 
scattering length changes sign in just one of the condensates. The coherent coupling is 
achieved by the quantum tunnelling through the central barrier of the double-well trap. 

There are two general geometric limits for the condensates trapped in a double-well 
potential: the one- and two-dimensional cases. In the one-dimensional case the condensates 
assume the cigar shaped form, while in the two-dimensional case they have the form of a 
pancake. We consider the two-dimensional case and rely on the important fact that the 
quasi two-dimensional Feshbach resonance is possible |l2|. The one-dimensional nonlinear 
Schrodinger equation describing the cigar shaped condensate has no collapsing solutions, 
they appear only when the transverse degrees of freedom are taken into account. 

The stationary states in a pair of tunnel-coupled BECs bifurcate from the zero solution. 
In this limit the two condensates with opposite interactions can be described by a single 
effective equation of the nonlinear Schrodinger type with the nonlinearity being repulsive 
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or attractive depending on the parameters. When the effective interaction is repulsive the 
system is stable, whereas, under certain condition on the number of atoms, the attractive 
effective interaction causes collapse. 

The surprising fact is, however, that the ground state in a pair of tunnel-coupled repulsive 
and attractive condensates can have a lower energy than the stationary state unstable with 
respect to collapse. Therefore, for some values of the parameters, there is an energy barrier 
for collapse in the system. Recall that in a single nonlinear Schrodinger equation it is, in 
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fact, the ground state which suffers from the collapse instability |9L llOl lll|. 

This work is related to our previous study ^18j of the tunnel-coupled condensates in two 
spatial dimensions, where we have considered the stationary states of the system but without 
analysis of the ground state. Here we focus on the analysis of the ground state in the system 
and its energy, whereas the detailed derivation of the coupled-mode system and some of the 
auxiliary analytical results can be found in our previous publications Refs. ^,19, 2ol |. 

The paper is organized as follows. Section O contains analysis of the ground states in the 
system in the coupled-mode approximation, i.e. for far separated wells of the double-well 
trap. We conclude with the summary of the results in section IIIII The details of the linear 
stability analysis of the stationary solutions are relegated to the Appendix. 



II. STATIONARY STATES IN THE COUPLED-MODE APPROACH 

Experimental realization of the scattering length management in just one of the two 
condensates trapped in a double-well potential necessitates that the wells be far separated. 
The lower two energy levels in such a double-well potential are quasi-degenerate and there 
is a large energy gap separating them from the rest of the energy spectrum. An ideal 
quantum gas in such a trap would occupy the degenerate subspace. For the condensate of 
an interacting gas there is a localized basis in the degenerate subspace which allows one to 
reduce the three-dimensional Gross-Pitaevskii equation to a system of linearly coupled two- 
dimensional equations (see also Ref. Indeed, due to the localization of the basis wave 
functions in the different wells of the double-well trap, the nonlinear cross-terms describing 
the atomic interaction between the two condensates are negligible. 

The central barrier in the double-well trap can be created by a laser modulation of 
the parabolic potential. Assuming that the parabolic trap allows for the pancake shaped 
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condensates, i.e. that the trap frequencies satisfy 7 = uJz/uj±_ 1, we can separate the 
spatial degrees of freedom, that is the order parameter for each of the two condensates 
factorizes: 2;, t) = ^^{x,y,t)ip\\{z), with tpw^z) being one of the wave functions of 

the locahzed basis. By the above factorization we neglect the nonlinear term in the Gross- 
Pitaevskii equation as compared to the kinetic term along the double-well trap, i.e. we 
assume 

where ^z is an estimate of the condensate length in the 2;-direction and g = Airas/m is the 
nonlinear coefficient in the Gross-Pitaevskii equation 1]. The longitudinal kinetic term is 



compensated by the external potential, hence iz = ^/h/ {rriLUz). 

The Hamiltonian of an atom in the double- well trap can be put in the form [1^ |Ji 

H = S\R){R\-Ki\R){L\ + \L){R\), (2) 

where \L) and \R) are the localized basis, S can be identified with the zero-point energy 
difference and K with the tunnelling coefficient. Under the condition supplemented by 
the assumption of far separated wells of the trap and the quasi degeneracy of the energy 
levels, i.e. K <^ fkJz and S ^ JtlOz, the Gross-Pitaevskii equation reduces to the coupled- 
mode system describing evolution of the condensates in the transverse dimensions: 

idru = (— + p^) u — \u\'^u — kv, (3a) 
IOtv = [-V"^ + p"^) V + {e + a\v\'^)v - nu (3b) 

(see also Ref. ^^). We have used the dimensionless variables: T = uJxt/2 and p = rjl^ 
(the coordinate in the transverse plane), with = ^Jfij (muj±). Here k, = 2K/(fiuj±), 
6 = 2S/{huj±), and a = ai^"* / 1 '^s"'' | is the relative scattering coefficient (we assume that 
ai"^ < and ai^'' > 0). The term in the system (jH)) is due to the common transverse 
parabolic trap. Note that due to 7 ^ 1 the coefficients k and e are not small in general. 
The transverse order parameters for the two condensates can be expressed as = 



{V^/i^)u and $^ = {VA/i^)v, with A = (87r|a^^|/ d^l^ii^-^ |l8|. For the current 
experimental traps the parameter A is on the order of 10^ — 10'^. In the numerics below we 
use the quantities Ny_ = J d^p |mP and = J d^p |f p referring to them as the number of 
atoms, for short, though the actual number of atoms is larger by the factor A. Finally, the 
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condition (^J in the dimensionless form reads: 

< -fRl aN, < -fRl, (4) 

where Ru and R^ are the radii of the two condensates measured in the units of i±. 

The stationary states with the chemical potential fi are sought in the form u = e~*^"^f/(p) 
and V = e~^^'^V{p). In general, there are two types of the real solutions U and V, satisfying 
the inequalities UV > (positive) and UV < (non-positive). The positive stationary states 
bifurcate from zero at fi = yUbif, where yUbif = 2 + e/2 — ^/e'^/A + k^, while the non-positive 
ones at /i[,j£ = 2 + e/2 + ^/e'^/A + k^. The bifurcation points are easily found by analysis of 
the linear part of the stationary coupled-mode system which follows from the system Q for 
the vanishing amplitude solutions U and V, given as U = Ae-p'l'^ and V = Be~p''^ with 
B = {2 — fi)A/K, and A —>■ 0. Being interested in the ground state of the system, we need 
to consider only the positive solutions since they bifurcate from the zero solution at a lower 
energy. Below the stationary solution (state) means the positive solution. 

The stationary states were found by using the numerical approach of Ref. [2^, i.e. 
by iterative solution in the polar coordinates of the nonlinear eigenvalue problem arising 
from the system Q (with fi being the nonlinear eigenvalue). Our method allows for the 
simultaneous computation of the stationary state and the chemical potential /i, with the 
same accuracy. We have used the Fourier pseudo-spectral method with the spatial grid 
containing 256 grid points. 

We have found that the sufficient condition for stability of the stationary state of the 
system © is given by the Vakhitov-Kolokolov (VK) criterion: dN/dfi < 0, with = N^+Ny 
being the total number of atoms. Moreover, if the effective interaction is repulsive (see 
below) then the ground state solution is unconditionally stable. The applicability of the VK 
criterion for stability in the case of the two-dimensional coupled-mode system is partially 
based on the numerical analysis of the appropriate eigenvalue problems, as discussed in the 
Appendix. 

If the tunnel-coupling were absent (k = 0) the attractive -u-condensate would collapse 
for Nu > Nth ~ 11.69. This threshold for collapse in the parabolic trap is g iven by the 
number of atoms in the so-called Townes soliton (see, for instance, Ref. |2l|). We have 
found that in the system Q with a > the stationary states unstable with respect to 
collapse always appear in the hmit of large negative and have — > Nth as /i ^ —00. For 
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some values of the system parameters the collapsing states exist for finite negative values of 
fi (see Fig. Q). As fi —oo they approach the Townes soliton in the u-component, while 
the f -component vanishes jl^. The numerically computed number of atoms as function of 
the chemical potential is illustrated in Fig. ^ For these values of the parameters k and a, 
the behavior of the number of atoms for positive values of the zero point energy difference is 
similar to that of figure ^a). It is seen that the curve = N{fi) develops a local maximum 
close to the bifurcation point as the zero-point energy difference decreases. This is related 
to the existence of a threshold value of the energy difference, given as Eq = h? /2 — 2, see 
below. 
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FIG. 1: The total number of atoms (decreased by the factor A) vs. the chemical potential in 
the units of fiw±/2 (the solid line) corresponding to a = 0.01 and k = 1. The dashed and dotted 
lines give the number of atoms in the u- and u-condensates, respectively. In panel (a) the unstable 
stationary state is attained in the limit of /i ^ — oo, while in panels (b) and (c) the unstable states 
appear for finite negative where dN/d^ > 0. The stable states correspond to dN/d^x < 0. 

Let us determine the energy of the stationary states. In the coupled-mode approximation 
the energy (measured in the units of Ahuj±/2) reads 



E 



dV <^ |Vm|^ + + p^(|m|^ + It^n + e\v\^ - k{uv* + vu* 
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FIG. 2: The energy of the stationary state (in the units of Ahuj±/2) vs. the total number of 
atoms (decreased by the factor A) corresponding to figure 1. The almost straight line in panel (c) 
contains the ground state of the system stable with respect to collapse. 

Using the coupled- mode system Q to express the kinetic energy of a stationary state one 
can reduce the expression for the total energy to the following 

E = ^N+^Jd^p{U^~ aV^} . (6) 

The numerically computed energy of the stationary state vs. the total number of atoms, 
corresponding to Fig. [TJ is illustrated in Fig. |H (to produce this figure we have computed 
the number of atoms and the energy separately by using the numerically found stationary 
states). It is seen that the collapsing states with N approaching the limit value A'th have 
the energy approaching zero. Note the change of sign of the derivative dE/dN at = as 
we decrease the zero-point energy difference in panels (a) through (c). (The fact that the 
zero-point energy difference is negative in figures |2fb) and (c) is a mere coincidence, see also 
figure m) 

The derivative dE/dN at N = can be found analytically as follows. For e = /ibif— /i 
the balance of the terms and fiU and the linear coupling of the U and V components 
in equation (jHaj) give U,V ^ VT^? while the characteristic radius of the solution is fixed by 
the trap size. Hence, from equation © we conclude that E = fiufN + 0{e^) due to the fact 
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that ~ |e| in this hmit. Thus OE/dN = fi^f + 0(e). 

If the sign of dE/dN at = is negative, that is e < £0 with 60 = /2 — 2, and 
the maximal value A^m of the number of atoms in a stationary state exceeds A^th then the 
ground state of the system with A^ > A^th is secured from collapse by an energy barrier. 
Indeed, first of all as either ji —00 or — » /ibif we have E ^ Hence, for /ibif < 
the function E = E{N) has a local minimum at A'm, where it assumes a negative value. 
Consider the stable stationary solutions with negative energy and the number of atoms 
A^th < N < (such solutions lie on the almost straight line of the curve E = E{N) in Fig. 
El^c)). The system in such a state can collapse if some finite amount of energy is passed to 
the condensates by an external perturbation, i.e. there is an energy barrier for collapse. 

In an attempt to understand the energy barrier for collapse let us first try to determine 
the effective interaction in the coupled-mode system. The effective interaction depends on 
the number of atoms and the double-well trap parameters. An example of the effective 
interaction switch from repulsive to attractive is given in Fig. |H1 which is a part of Fig. ^c) 
zoomed in about /ibif- The interaction changes sign at dfi/dN = when a zero eigenvalue 
appears the linear stability spectrum (see Eqs. ()A2|) - ()A4|1 in the Appendix; the positivity of 
the linear stability spectrum is a sign of the effective repulsion in the system, see below) . As 
shown in the Appendix, the stationary states are unconditionally stable when the effective 
interaction is repulsive. 




-1.303 -1.3025 -1.302 



FIG. 3: Zoomed in part of figureQ^c) about /i = /ibif which manifests the effective repulsion in the 
coupled- mode system in the vicinity of the bifurcation point. Here the axis units are as in figure ^ 

The effective interaction can be analytically studied as e = /ibif — /i — >■ 0. First we use 
equation ()3b|) to express V in terms of U and its derivatives as follows 

V = K(3U - an^P^U^ + Kp^{V^ + 2 - p^)U + 3a\^(3'^U^ + ..., (7) 
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where (3 = (e + 2 — and the higher order terms are represented by dots. Equation 
(|7j) is derived in a similar way as an analogous equation for the one-dimensional coupled 
mode system in Ref. 19|, one only has to replace the second-order differential operator by 



+ 2 — and p by p — 2. Following Ref. jiol we substitute the expression ((7j) into equation 



aj) and get 



where /xq = 2 + e/2 + ^/eyAT~i^. 

In the limit of the small amplitude solution the interaction sign is defined by the coefficient 
at the cubic term (apart from the exceptional case of zero). The effective interaction is 
attractive for e > Sint with eint = K(a^^^ — a~^/^). This formula accounts for Fig. 01 since 
there we have e < ei^t = —2.85. 

Comparison of Figs.|2fc) and|Slleads to the conclusion that the energy barrier for collapse 
in the system cannot be explained by the repulsive effective interaction. Indeed, the interval 
of the axis occupied by the stable ground states in Fig. |2fc) is much wider that of the 
repulsive interaction in Fig. El 

Moreover, we get ^int < £o for k, > Kcr, where = a/Q^ + 4 — Q with Q = a~"^/^ — a"^/^. 
Hence, for k > and eint < s < Eq there are stable ground states with negative energy in 
the limit p /ibif, while the effective interaction is attractive. An example is provided in 
Fig. m where e = 1 while ^int = —10.65 and = 5. 
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FIG. 4: An example of the energy barrier for collapse in the system with the attractive effective 
interaction. In this case Sint = —10.56, while £q = 5. Here the axis units are as in figureEl 

Fig. m also is a counterexample to the possible explanation of stability of the ground 
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state with a negative energy by the fact that the repulsive condensate has a lower zero-point 
energy {e < 0), i.e. is placed in the lower well of the double- well trap. 

To understand the existence and stability of the stationary states with negative energy let 
us examine the expression for the energy (jU)). As shown above, for small number of atoms, 
i.e. for /ibif, the linear part of the r.h.s. of equation (jU)) is dominant. But for a large 
number of atoms the nonlinear part must dominate. In this case, it is seen that the energy 
decreases with increase of the number of atoms if the term aV^ is dominant in the integral 
(it is important to note that the two terms in the integral on the r.h.s. of formula (jH)) are due 
to an interplay of the atomic interactions, the kinetic energy, the quantum tunneling and 
the zero-point energy difference in the system, since we have expressed the kinetic energy 
of the condensates in terms of the other energies to obtain formula ©). Indeed, comparing 
the figures H and 121 we see that when dE/dN < a much larger share of atoms is confined 
in the repulsive condensate (figures HJ^c) and|2tc)). Therefore, such a stationary state does 
not collapse, since the energy would increase if some of the atoms tunnel from the repulsive 
condensate to the attractive one. This causes the energy barrier for collapse. 

The stable ground states with negative energy in the system of tunnel-coupled repulsive 
and attractive two-dimensional condensates are analogous to the "unusual solitons" in the 
o.e-d,.en.o.al case Q wMch also appea. Cose to the bifu.eation po,„t. The unusual 
soliton is a bright soliton which has a much larger share of atoms confined in the repulsive 
condensate. 

III. CONCLUSION 

In conclusion, we have found the ground states in the system of tunnel-coupled repulsive 
and attractive condensates, which are stable with respect to collapse. The coupled-mode 
approximation have been used and is justified by the condensates being far separated in the 
wells of the double-well trap. We predict the existence of an energy barrier for collapse in 
the system, which is related to the appearance of the stable ground state with a negative 
energy. The condensates in such a ground state can collapse only if they are externally 
excited with a finite amount of energy being passed to the system. 

The ground state with a negative energy is characterized by a large share of atoms trapped 
in the repulsive condensate. The contribution to the energy from the attractive condensate 



10 



is positive, whereas from the repulsive one is negative (the contribution of each condensate 
is determined as the resuh of an interplay between the kinetic energy, the interaction energy 
and the zero-point energy, additionally there is the coupling energy due to the quantum 
tunneling). Therefore, the energy increases if some atoms tunnel from the repulsive to 
the attractive condensate. This explains the existence of the energy barrier for collapse in 
the system. The energy barrier appears for a wide region of the parameters space under 
the condition that the zero-point energy difference is less than a threshold value. The 
latter depends on the tunneling coefficient and can be both positive and negative (i.e. the 
condensate with repulsion can be placed in the upper well of the asymmetric double-well 
trap). 

As is known, the Gross-Pitaevskii equation is sufficient to describe stable stationary 
states of the Bose-Einstein condensates at zero temperature, due to the similarity between 
the Bogoliubov-de Gennes equations and the equations describing evolution of a linear per- 
turbation of the condensate order parameter 22|. Therefore, at zero temperature, given 
that the cloud of non-condensed atoms can be neglected, the energy barrier for collapse 
phenomenon can be observed. 
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APPENDIX A: LINEAR STABILITY ANALYSIS 

We will consider only the axially symmetric stationary points of system (jS)): u = 
e~^^^U{p) and v = e~'^^'^V{p). Stability can be established by considering the eigenvalue 
problem associated with the linearized system. Writing the perturbed solution as follows 

u = e-'^^ {U{p) + e-'^^W(p)} , V = e-'^^ {V{p) + e-*^^V(p)} , (Al) 

where (W,V) is a small perturbation mode with frequency fi, one arrives at the following 
linear problem for the eigenfrequency: 




with 

Here the scalar operators are defined as follows 



lS") = -(V2-3aV^2^p-£)+p2. (A4) 

First of all, the matrix operator Aq is non-negative for positive stationary solutions, i.e. 
satisfying UV > 0. Indeed, the scalar operators on the main diagonal of Aq can be cast as 
follows 

what can be easily verified by direct calculation. Therefore the scalar product of Aq with 
any vector X = (Xi(p), X2(p)) is non- negative: 

(X|Ao|X) = J dV |x*4")Xi + X*4")X2 - fi:(X*X2 + X*Xi)} 



>K / dV 



.2 



> 0. 



Here we have used the positivity of the operators — ^Vf/^V^ and —yVV'^Vy. The operator 
Aq has one zero mode given by the stationary point itself: Z = {U,V), AqZ = 0. Non- 
negativity of Aq is an essential fact for the following. Thus the non-positive solutions, i.e. 
satisfying the inequality UV < 0, are discarded from the consideration. 

The lowest eigenfrequency of the linear stability problem can be found also by minimizing 
the following quotient 

. r.in^^j^n. (A5) 

in the space orthogonal to the zero mode of Aq: {Z\X) = (here {X\Y) = J d^p(X*Y'i + 
12-^2))- Equation ()A5|) follows from the eigenvalue problem rewritten as AqAiX = Q'^X 
with X = {Ur,Vr). 

The imaginary eigenfrequencies fl, which mean instability, appear due to negative eigen- 
values of the operator Ai. If there is just one negative eigenvalue, then the Vakhitov- 
Kolokolov (VK) stabihty criterion dN/dfi < applies, which can be established by a simple 
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repetition of the arguments presented, for instance, in Ref. [20|. The hmit on the number of 
negative eigenvalues is related to the fact that the minimization of the quotient in equation 
()A5|) is subject to only one orthogonality condition, thus only one negative direction can be 
eliminated by satisfying this orthogonality condition. 

To apply the VK criterion for stability of the stationary points in the 2D case (as compared 
to the ID case j^]) one has to rely on numerics to establish the number of negative eigen- 
values of the operator Ai. The following simple strategy was used. The eigenvalue problem 
was reformulated in the polar coordinates (p, 6) and the operators were expanded in Fourier 
series with respect to 9 by the substitution: ^ Vp — n^/p^, where Vp = dp + p~^dp. 
Noticing that the orbital operators Ai^„ = Ai(V^ — > — / p^) are ordered as follows 
^i,n+i > Ai „, we have checked for the negative eigenvalues of the first two orbital operators 
with = 0, 1. It turns out that Ai^i is always positive, while the operator Ai^o has one 
negative eigenvalue or none (the latter corresponds to the effective repulsive interaction in 
the system). Thus the VK criterion applies in the former case, while in the latter one the 
stationary point is unconditionally stable. 
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